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ial An abbreviated version of this protocol was published in eLIFE in Sep 2019 
Molecular function limits divergent protein evolution on planetary timescales 


DOI: 10.7554/eLife.39705 


Detailed protocol 


Inferring the divergence rate between orthologs as a function of divergence time or protein distance 


1. Define tand y vectors to fit divergence model. For each pair of orthologs, values in t are the divergence time between species expressed in billions of 
years, values in y are the percent sequence identity between corresponding orthologs. 


2. Use the fitnim function in matlab (or an alternative nonlinear regression method) to fit equation 3 from Konate eft al. and retrieve the values for the Rg and 


alpha parameters. (relevant matlab code shown below). 
* a 
y= 100*(#2* +1) 


model3=@(b,x)100.*((b(1).*x/b(2))+1).*-b(2) 
beta=[1 0.5]; 

fit=fitnim(t,y, model3,beta); 

[RO,alpha]=fit.C oefficients.Estimate; 


3. Use the obtained Rg and alpha parameters in equation 5 from Konate et al. to determine the divergence rate (dy/dt) at an arbitrary t expressed in billions 


of years. 
dy Ro*t 
o _ 100 *R (#2 + 1) 


4. Alternatively, use the obtained Rg and alpha parameters in equation 4 from Konate et a/. to determine the divergence rate an arbitrary protein distance 


expressed as a fraction (i.e. 1 — percent identity/100) 


SP = Ro%(1— dD)" 
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